Density correlations in the half-filled Hubbard model 
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We consider density-density correlations in the one-dimensional Hubbard model at half filling. 
On intuitive grounds one might expect them to exhibit an exponential decay. However, as has 
been noted recently, this is not obvious from the Bethe Ansatz/conformal field theory (BA/CFT) 
approach. We show that by supplementing the BA/CFT analysis with simple symmetry arguments 
one can easily prove that correlations of the lattice density operators decay exponentially. 



Recently, density correlation functions have been stud- 
ied for various multicomponcnt one-dimensional elec- 
tron systems within the framework of the bosonization 
approactoa. An important issue raised in these works 
is the question whether or not these correlators exhibit 
a power-law behaviour at large distances in the presence 
of a Mott-Hubbard gap (see also Refs. H||). For the par- 
ticular case of the half- filled, one-dimensional Hubbard 
model a bosonization analysis yields the expected expo- 
nential decay, which at first sight appears to be at-,odds 
with the results obtained from the Bethe AnsatzB (see 
e.g. [20] of Ref. |). 

In one spatial dimension, the asymptotic behaviour of 
correlation functions can be analyzed by means of the 
BA/CFT approach. Furthermore, all possible critical di- 
mensions can be obtained from the exact solutionB. Since 
the expansion of the lattice operators in terms of the 
conformal fields is not known explicitly, the resulting ex- 
pressions contain unknown amplitudes. In particular, the 
leading term in the BA/CFT expression for the correla- 
tion function under question could actually vanish. This 
is the origin of the controversy mentioned above. In cer- 
tain cases, including the one discussed here, it is possible 
to resolve this issue by employing symmetry considera- 
tions of the underlying microscopic model. 

In this note, we are interested in connected correlation 
functions of electron densities at half filling. Using the se- 
lection rules of conservation of numbers of electrons with 
spin up and down, one can show from the exact solution 
(specialized to the case of zero magnetic field) that 
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where v is the spin velocity, ao is the lattice spacing and 
Ai t 2 are unspecified amplitudes that depend on the phys- 
ical parameters of the problem i.e. interaction strength 
and in the general case magnetic field. Note that in (|l|) 
we have omitted similar, subdominant terms. 

We will now show that by using the full set of selection 
rules for the half-filled Hubbard model, one can establish 
that the amplitudes of all contributions to (0), that decay 



algebraically, are identically zero and that the correlator 
thus decays exponentially. 

The Hamiltonian of the half-filled Hubbard model is of 
the form 



(2) 



Here c^, cj are canonical annihilation and creation op- 
erators for electrons with spin a on site i, rii = — 

J2 a c L c i<>~ i s t ne number operator density for electrons. 

In addition to the SU (2) invariance in the spin degrees 
of freedom the Hubbard Hamiltonian on a bi-partitc lat- 
tice commutes with the generators, of a second SU(2) 
related to particle-hole symmetryQij. In the one dimen- 
sional case the local generators of this ^-pairing SU (2) 
are 



ifCj-lj Vj : 



(4). 
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satisfying commutation relations [rjfiVj] = an d 
[Tjj,TjJ'] = 2?7|. Combining the two symmetries one 
obtains the well-kown SU(2) ® SU{2)/Z 2 = SO{A)- 
symmetry of the Hubbard model. 

At half filling the ground state |0) is an SO(4)- 
singletEm. Denoting by n(t, x) the density operator at 
x = a oj and time t we want to determine the asymp- 
totic behaviour for j — » oo of the connected correlation 
function 

G nn (t,x) =4(7f(t,:r)7f (0,0)) 

= 4^(f7|,7 z (t,x)|n)H^(0,0)|fi) 
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Hence the matrix elements 
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determine which intermediate states will contribute to 
the asymptotic behaviour of (Q). 

In the BA/CFT approach to correlation functions 
in integrable models one usually proceeds as follows: 
Firstly, one calculates the finite-size scaling behaviour of 
the energies e n and momenta k n of the low-lying states 
entering the expression (Q) . Secondly, one uses the quan- 
tum numbers of the operator entering the matrix element 
to eliminate "incompatible" intermediate states. Finally, 
the CFT relation between finite-size energies and scaling 
dimensions is used to obtain the asymptotic (power law) 
behaviour of correlation functions by summing over the 
remaining low- lying intermediate states. 

In Ref. | only particle numbers for spin-cr electrons 
were used in the second step. This corresponds to tak- 
ing into account only quantum numbers associated with 
the Cartan generators S z and rf but not those with the 
total spins S 2 and fj 2 . As a result one obtains (|l|), which 
seems to indicate that G nn decays algebraically for large 
distances in spite of the presence of a charge gap at half 
filling. We will now show, by considering the full set of 
SO (A) quantum numbers, that only intermediate states 
with a gap contribute to (j^) and the density correlations 
decay exponentially. 

Let us calculate the total ry-spin quantum number of 
the state T]* using the commutation relations between 
the ry-pairing operators (and the fact that \ is a singlet 
under the ^-pairing SU(2)) one easily obtains 
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showing that this state is a triplet of the ^-pairing SU(2). 
As a consequence, only intermediate states which are r]- 
pairing triplets can contribute to the correlation function 

As is well known only pure spin excitations are gap- 
less in the repulsive half filled Hubbard model. However, 
these are all singlets of the ^-pairing ST/ (2)0. Therefore 
the corresponding matrix elements (ra) vanish identically 
and cannot contribute in G nn . We conclude that the 
lowest energy intermediate states with non-zero matrix 
elements (||) are holon-antiholon scattering states with 
energy above the Mott-Hubbard gap. As a result the 
density-density correlation function exhibits exponential 
decay at large distances for any positive U 



G nn {t,x) — ► exp(-ax) , 



oo. 
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The determination of a(U) is an interesting open prob- 
lem. On general grounds we expect4t-to be proportional 
to the holon gap, which is given byl!30 (see also Ref. H2) 



A(*7) = --2 + 4^(-l) n y/1 + (nU/2) 2 - nU/2 



VUexp{-2ir/U) , for U -> . (8) 

If one considers correlation functions of "point-split" 

densities such as = J^a- c ]<7 C j+i,v tne above sym- 

metry argument does not imply the vanishing of e.g. the 
matrix element (ST|Mp s )(0, 0)|0), where |ST) denotes a 
spin-triplet excitationtS. Given that symmetry does not 
force such matrix elements to vanish we expect them to 
be nonzero, which then immediately implies an algebraic 
decay like in (|l]). 

We emphasize, that to establish the exponential decay 
of G nn in the above way it is essential that: (1) there 
is an exact symmetry in the charge sector of the micro- 
scopic Hamiltonian (j^). In general this symmetry may 
not be obvious (as is the case for the Hubbard model). 

(2) the ground state is a singlet of the corresponding alge- 
bra. Note that the above considerations still hold in the 
presence of a magnetic field as the ground state (in the 
half-filled band) remains a singlet of the ^-pairing SU(2). 

(3) All charged (non-singlet) excitations are gapped. 
We note that these conditions are fulfilled for the half- 
filled Hubbard model on a bipartite lattice in any dimen- 
sion, providedpthat U is larger than the critical Mott- 
Hubbard valueS. 

Analogous conclusions can be reached for spin-spin cor- 
relations in the attractive Hubbard model by employing 
discrete symmetries of the Hamiltonian (|^). Under the 
particle-hole transformation for spin up 
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the Hamiltonian transforms according to H(U) — > 
7i(—U), whereas ^-pairing and spin SU(2) symmetries 
are interchanged S a <-*■ rf, a = +, — , z. Furthermore the 
ground state of the attractive Hubbard model is a spin 
singlet. This implies that spin-spin correlation functions 
in the attractive Hubbard model decay exponentially at 
large distances 



(S z (t, x)S z (0, 0)) — ► exp(-/3i) , x 



(10) 



where (3 > for any U < 0. 
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